Pattern formation with a conservation law 



P. C. Matthews and S. M. Cox 

School of Mathematical Sciences, University of Nottingham, 
University Park, Nottingham NG7 2RD, UK 

Abstract. Pattern formation in systems with a conserved quantity is considered 
by studying the appropriate amplitude equations. The conservation law leads 
to a large-scale neutral mode that must be included in the asymptotic analysis 
for pattern formation near onset. Near a stationary bifurcation, the usual 
Ginzburg-Landau equation for the amplitude of the pattern is then coupled to 
an equation for the large-scale mode. These amplitude equations show that for 
certain parameters all roll-type solutions are unstable. This new instability differs 
from the Eckhaus instability in that it is amplitude-driven and is supercritical. 
Beyond the stability boundary, there exist stable stationary solutions in the form 
of strongly modulated patterns. The envelope of these modulations is calculated 
in terms of Jacobi elliptic functions and, away from the onset of modulation, is 
closely approximated by a sech profile. Numerical simulations indicate that as 
the modulation becomes more pronounced, the envelope broadens. A number of 
applications are considered, including convection with fixed-flux boundaries and 
convection in a magnetic field, resulting in new instabilities for these systems. 



PACS numbers: 47.54.+r, 02.30.Jr, 05.45.-a 

1. Introduction 

The mathematical theory of pattern formation has a wide range of applications. In 
the field of fluid mechanics, the most widely studied example is Rayleigh-Benard 
convection and its variants, but there are many other examples such as Taylor-Couette 
flow and solidification of alloys. Pattern-forming systems can also be found in many 
chemical and biological systems (see the review by Cross and Hohenberg 1993). 

In this paper we consider pattern formation in one space dimension in systems 
sufficiently large to allow a large-scale modulation of the pattern. It is assumed that 
the system has some reflection symmetry, that the onset of instability is stationary, and 
that the onset of pattern formation occurs with a non-zero wavenumber. In this case, 
the equation governing the modulation of the small-amplitude pattern is in general 
the real Ginzburg-Landau equation. This well known result has recently been given 
rigorously by Melbourne (1999). The behaviour of the real Ginzburg-Landau equation 
is straightforward: a regular periodic pattern is stable provided its wavenumber lies 
within the Eckhaus band (Eckhaus 1965); patterns outside this band are unstable, and 
the instability is subcritical (Fauve 1987, Tuckerman and Barkley 1990). Physically, 
patterns in which the wavenumber is either too large or too small adjust to a different 
wavenumber by adding or losing waves. 

However, there is a large class of problems in which the onset of pattern formation 
at finite wavenumber is not governed by the Ginzburg-Landau equation. Such 
problems arise when the system possesses a conservation law (Cross and Hohenberg 



Pattern formation with a conservation law 



2 



1993, p. 883). To be specific, we suppose that this conservation law may be written 
in the form 



with appropriate boundary conditions, from which it follows that w is conserved, that 
is, 



where the integral is taken over the entire spatial domain. The conservation law leads 
to a neutral large-scale mode, which must be included in the dynamics if a correct 
description of the behaviour near onset is to be obtained. 

The resulting behaviour depends on how the reflection symmetry acts on the 
large-scale mode. The case in which the large-scale mode changes sign upon reflection 
arises when the system has Galilean invariance (Coullet and Fauve 1985). This case 
was studied by Matthews and Cox (1999), who showed that all steady patterns are 
unstable, with the instability growing more rapidly than the rate of formation of the 
pattern. As a consequence, a highly disordered state is found, even near onset. 

In this paper we study the case where the large-scale mode does not change sign 
under reflection. Physically, this corresponds to the mode representing a density rather 
than a flow. Similar problems have been considered by Riecke (1992a, 1992b, 1996) 
and Herrero and Riecke (1995) for convection in binary mixtures, but in this case the 
bifurcation is oscillatory and asymptotically self-consistent amplitude equations were 
not obtained. 

Experiments demonstrating the effects of coupling between short-scale pattern 
modes and large-scale mean flows have been carried out in a variety of flows. In 
surface tension-driven convection, VanHook et al. (1995, 1997) have shown how the 
presence of a large-scale surface-deformational mode can lead to dry spots and high 
spots. In a two-layer viscous shear flow (Barthelet and Charru 1998, Charru and 
Barthelet 1999), finite-amplitude interfacial waves become unstable through coupling 
to a large-scale mode which may be traced back to conservation of mass of each 
fluid. In liquid crystals (Hidaka et al. 1997), complicated spatio-temporal convection 
patterns are traced to the presence of a large-scale mode arising in the degeneracy of 
the system. 

The paper is organised as follows. In § |^ we derive the governing amplitude 
equations by symmetry and scaling arguments and then, in § || analyse their 
properties. We consider the stability of small-amplitude solutions, and the behaviour 
of fully nonlinear solutions of the amplitude equations, in particular localised states. 
In § ^ we illustrate how the insights we have gained apply to a model partial 
differential equation (PDE), through numerical simulation. We also demonstrate 
how the behaviour of the full PDE diverges from that predicted by the leading-order 
amplitude equations, as the parameters are varied, and analyse this divergence in § ^|. 
We then describe, in § ^], three physical examples in which our amplitude equations 
arise. In § |t] we summarise our results. 

2. Derivation of mean-mode and pattern-mode equations 



dw d „. , 




w dx = 



In this section we derive the governing amplitude equat ions for pattern formation with 
a conservation law, using two different methods. In § 2T the equations are deduced 
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from symmetry considerations, while in § 2.2 a model PDE is introduced and the 
amplitude equations are derived by asymptotic expansions near the onset of pattern 
formation. 

2.1. Derivation of equations from symmetry 

Consider a PDE or system of PDEs with a conserved quantity and a basic zero state. 
If a quantity w(x,t) is conserved, it must be possible to write dw/dt as the divergence 
of a flux. It follows that when the expression for dw/dt is linearised, there can be 
no terms linear in w except those involving cc-derivatives. If we assume, furthermore, 
that the system possesses the reflection symmetry 

x — > — x, w — > w (2-1) 

then only even ^-derivatives occur. Thus w is density-like (a velocity-like quantity 
would satisfy x — > —x, w — > — w instead; see Matthews and Cox (1999), Tribelsky 
and Velarde (1996)). On large spatial scales, the leading linear term is the second 
derivative w xx , so w{x, t) obeys the diffusion equation and the growth rate A of modes 
with wavenumber k is proportional to k 2 . This indicates that large-scale modes are 
almost neutral, evolving on a long timescale, and must be included in the analysis of 
pattern formation. 

Suppose that the system supports stationary patterns with a wavenumber k, 
which we take to be the critical wavenumber according to linear stability theory; we 
may set k = 1 with no loss of generality. Near the onset of pattern formation the 
appropriate ansatz is then 

w{x,t) ~ A(X,T)e ix + A*(X,T)c- ix + B{X,T), (2.2) 

where X and T are rescaled forms of x and t, and A and B are small amplitudes. 
Note that the amplitude A is complex, with A* denoting the complex conjugate of A, 
but B is real. 

The amplitude equations for A and B can be written down by imposing the 
symmetries of translation in x 

A^Ae lx , B^B (2.3) 

and reflection 

X-y-X, I -*•!*, B^B. (2.4) 

In the absence of the large-scale mode B, A obeys the real Ginzburg-Landau equation. 
Coupling terms in this equation must include A because of (^j), so the leading 
coupling term in the A equation that is consistent with the symmetries is AB. The 
coefficient of this term is forced to be real by ( 2.4). I n the absence of coupling terms, 
the equation for B is the diffusion equation, and (2.3) forces coupling terms to involve 
equal numbers of A and A* terms. Since B is conserved, all terms in the B equation 
involve derivatives. The leading coupling term in the B equation that is consistent 
with (U) is therefore {\A\ 2 ) XX . 



These considerations give rise to amplitude equations of consistent asymptotic 
order if we introduce a small parameter e such that A = eA, B = e 2 B, X = ex, 
T = e 2 t. Furthermore, we take the system to be a distance of 0(e 2 ) from the marginal 
surface in parameter space. With these scalings, the amplitude equations truncated 
to leading order are then 

A T = A + A XX -A\A\ 2 -AB, (2.5) 

B T = aB xx +^\A\ 2 ) xx , (2.6) 
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Figure 1. The growth rate A of the linear modes of (2.7) as a function of 
wavenumber k, with r = 0.1. 



where we have assumed that, in the absence of B, the bifurcation of the pattern mode 
is supercritical. All coefficients in (2J5) have bee n set to unity by rescaling T, X, A and 
B. There are two coefficients, a and /i, in ( |2.6|) that cannot be removed by rescaling. 
We impose the condition a > to ensure that the long-wavelength modes are weakly 
damped in the absence of A, but /i may have either sign. Equations analogous to 
( [2.5| )-( |2~6| ) have previously been derived in the context of the secondary stability of 
cellular patterns (Coullet and Iooss 1990). 



2. 2. Derivation of equations from a model PDE 
Consider the PDE 



dw 
~di~ 



dx 2 



rw — I 1 



dx 2 



w — sw 2 — w 3 



(2.7) 



where the terms inside the brackets are just those of the much-studied Swift- 
Hohenberg equation (Swift and Hohenberg 1977) with a symmetry-breaking quadratic 
term (Sakaguchi and Brand 1996, Matthews 1998). This equation has the required 
properties of reflection symmetry and a conservation law for w(x, t). 
The growth rate A of modes of wavenumber k in (|]t]) is given by 

\ = k 2 {r-(l-k 2 ) 2 ), (2.8) 

so for < r <C 1 a narrow band of wavenumbers near k = k c = 1 is unstable (figure |l|) . 

We now follow the standard asymptotic methods to describe the behaviour of 
( |2.7| ) for small r, by writing 

r = e 2 r 2 , T = e 2 t, X = ex. (2.9) 

The appropriate asymptotic expansion for w(x,t) is then 

w(x, t) = eA(X, T)c ix + eA* (X, T) C -' 1X + e 2 B(X, T) 

+ e 2 C(X,T)e 2ix + e 2 C*{X, T) e - 2ix + 0(e 3 ), 



(2.10) 
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where the large-scale mode B has been introduced at order e 2 . 

The expansion (2.10) is substituted into (2.7) and the system is then solved at 
successive orders of e. At 0(e), ( pT7| ) is automatically satisfied. At 0(e 2 ), terms 
involving e 2lx appear, and from these we deduce that C — — sA 2 /9. At 0(e 3 ), the 
amplitude equation for A is deduced from the terms in e lx . The equation for B follows 
from the terms independent of x at 0(e 4 ). These equations are 

A T = r 2 A + AA XX - (3 - 2s 2 /9)\A\ 2 A - 2sAB, (2.11) 
B T = B xx + 2s(\A\ 2 ) xx . (2.12) 

The sign of s is thus immateri al, an d the bifur cation is supercritical if s 2 < 27/ 2. If 
this condition is satisfied th en ( 2.11 ) and ( 2.12 ) can be rescaled to give (2.5-2.6). 

Although the system ( |2.5| )-( |2.6"1 ) is generic for pattern-forming systems with a 
conservation law, in fact a conservation law alone is not sufficient to yield these 
equations. Only the first of the two nonlinear terms in (|]^) generates the coupling 
term in (2.5), so if we set s — there is no coupling and the dynamics is then governed 
by the Ginzburg-Landau equation alone. Furthermore, if the original PDE has the 
symmetry w — > — w then all nonlinear terms in (2.11) are cubic and again no coupling 



terms arise. This restriction will apply when we consider physical examples in 



3. Analysis of the amplitude equations 



In this section we analyse the steady solutions of the amplitude equations (2.5 
2.6). We begin by considering small-amplitude, periodic solutions (rolls). We then 
investigate the super/subcritical nature of the instability, and, beyond the point where 
rolls become unstable, we calculate fully nonlinear solutions of (|2.5| |2.0|). 



3.1. Stability of patterns 

In this section we analyse the stability of periodic solutions in the form of rolls. Rolls 
with wavenumber k c + eg correspond to solutions of (2.5-2~6|) of the form A — aoe lqX 



and B = 0, where the amplitude ao may be taken to be real and to satisfy a 2 , = 1 — q 
To analyse the stability of these rolls, we set 

A= {a Q + a(X,T))e i<!X , B = b(X,T), 

and linearise (2.5-^76 ) in a and b. Upon setting a — V(T)e llx 



b = [/(T)e 1 



IX 



U*(T)e- ilx , we find 
V = - (2q + 1)IV - a 2 (V + W)- a U, 
W = - {-2q + 1)IW - a 2 {V + W)- a U, 
U = -ol 2 U - tm Q l 2 {V + W). 



W*(T)e 



(3.1) 
ux and 

(3.2) 
(3.3) 
(3.4) 

Solutions of (3.2-3.4) proportional to e AT have a growth rate A that satisfies the cubic 
equation 

A 3 + (2l 2 + al 2 + 2al)\ 2 - I 2 (2^ag - 2aal - 2a 2 , - I 2 - 2al 2 + 4g 2 )A 

-Z 4 (2^a 2 - 2aa 2 ~ al 2 + Aaq 2 ) = 0. (3.5) 

For modes with A real, there is instability if 

2(n - a)a 2 + o-(Aq 2 - I 2 ) > 0. (3.6) 
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The most dangerous disturbance therefore arises in the limit I — > 0, and so the criterion 
for instability is 

- a)al + 2crq 2 > 0, (3.7) 

or, equivalently, 

Li(l-q 2 )>a(l~3q 2 ). (3.8) 



From (3.7) it is clear that all rolls are unstable if 

fx > a. (3.9) 

If (i < a then rolls of sufficiently small amplitude are unstable; the criterion for 
instability can be expressed in terms of q as 

q 2 > (<7-/i)/(3<r-/z). (3.10) 

It is clear that when the coupling between the mean and pattern modes is removed 
(/x = 0), or more generally when \cr/fx\ ^> 1, the usual Eckhaus criterion q 2 > 1/3 for 
instability is recovered. In general, the band of stable rolls is narrowed if fx > and 
widened if fx < 0. 



Hopf bifurcations are also permitted in (3J5). We are able to make analytical 
progress in the smalW limit, where one solution of (3.5) is A = — 2a 2 , while the other 
two obey 

\ 2 a 2 Q + M 2 {al(l + a - fx) - 2q 2 ) + Z 4 (cr -fx- 3q 2 a + q 2 fx) = 0. (3.11) 

It follows that there is oscillatory instability if q 2 > (1 + a — fx)/ (3 + a — fx) provided 
that fia — a 2 — fi > 0. From these conditions it can be deduced that a Hopf bifurcation 
is possible only if a < 1, fi < and q 2 > 1/3. There is a Takens-Bogdanov point (a 
repeated zero eigenvalue) at q 2 = 1/(3 — 2a) when fj,(a — 1) = a 2 . 

The stability res ults are summarised in figure |[ Note that the stationary 
bifurcation given by ( |3 .10 ) depends only on the ratio fi/a but this is not true for 
the Hopf bifurcation. 



The mechanism of instability in (2.5-2.6) is rather different from that of the 



Eckhaus instability. Even rolls with q = are unstable for fi > cr, and in this case 
the unstable mode has V — W so that a is real and the instability is amplitude- 
driven rather than phase-driven. If there arc regions where \A\ 2 has a minimum, then. 



according to (2.6), B grows in these regions if fx > 0. This leads to a positive value of 



B which then reduces the size of A further, through the coupling term in (|2 
3. 2. Nonlinear development of the instability 



The results of § 3.1 show that for fx > a, all periodic patterns in (2~5]-2~6|) are unstable 



We now consider the nonlinear development of the instability, to determine whether 
the bifurcation is subcritical or supercritical. The Eckhaus instability, for example, is 
known to give rise to a subcritical bifurcation (Fauve 1987, Tuckerman and Barkley 
1990). 

For simplicity we consider the case q = 0, corresponding to rolls with wavenumber 
k c . These are the last rolls to become unstable as ix is increased. In this cas e th e 
stationary solution is A = 1, B = 0. Writing A = 1 + a(X, T), B = b(X, T), cf. |y]), 
we find that the imaginary part of a decays so that we may take a to b e real. The 
nonlinear equations for the real- valued quantities a and b are then, from ( |2.5| |2.6| ), 

a T = -2a + a X x - 3a 2 - a 3 - b - ab, (3.12) 

b T = ab X x+fx(2a + a 2 ) xx . (3.13) 
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Figure 2. Region of stable rolls in ( |2,E| - [2,q ), The solid line shows the stationary 
bifurcation, which is valid for any <r, and the dotted line shows the Hopf bifurcation 
for cr = 1/4. For fi > a all rolls are unstable, while for fi < a there is a band of 
stable rolls that widens as /j, is decreased. The usual Eckhaus result of stability 
for q 2 < i corresponds to fi = 0. 



A weakly nonlinear analysis is now employed, by writing 
/i = /i + e 2 M2, T 2 = e 2 T, 

a = eat + e 2 a 2 + 0(e 3 ), b = eh + e 2 b 2 + 0(e 3 ), 
where the and bi are functions of X and T 2 . By considering a single mode, 
o-i = an (22) cos IX, b\ = buiTz) coslX, 



we find that the terms of order e in (3.12-3.13) give 



6n + (2 + l 2 ) ail = and fi = (1 + l 2 /2)a, 
so the instability condition (3.9) is recovered for small I. At 0(e 2 ) we deduce that 



«2 = «20 + «22 cos 21X, b 2 = b 22 cos 21X, 

where 

«20 = (i 2 - l)a?i/4, a 22 = fln/4 and 6 2 2 = -(1 + l 2 /2)a\ v 

Applying the solvability condition at 0(e 3 ) gives the equation for the evolution on the 
slow timescale T 2 : 

, ,9 „ ,9xdon al A (l 2 — 1) o 

(al 2 + 2 + Z 2 )— ^ = 2ZV 2 an + 4 -o? x . 

dJ 2 z 

Hence the bifurcation is a supercritical pitchfork provided that I < 1. In a large 
periodic domain, the first modes to become unstable have small /, and so the 
bifurcation will in practice be supercritical. 
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3.3. Nonuniform steady solutions 

We now extend the analysis of the previous subsection to consider fully nonlinear 
stationary solutions of the amplit ude equations ([T^)-( |2~6|) . Writing A = i?cxpiO, we 
find from the imaginary part of ( |2.5| ) that (R 2 &x)x = and so @x = fti? -2 , where 
ft is constant. From ( [2.6| ) we find, for solutions periodic in X, that 

B = n'((R 2 ) - R 2 ), (3.14) 



where (• • •) denotes a spatial average. In deriving (3.14), we have taken (B) = 0, 
which corresponds to an initial condition such that (w) = (and hence (w) = for 
all time). Any other initial value for (w) is equivalent to (w) — under a rescaling of 
the parameters r and s in fl2.7D - From the real part of ( |2.5| ), it then follows that R 
satisfies 

= (1 - n'{R 2 ))R + R X x - (1 - v')R 3 - ft 2 iT 3 . (3.15) 

This equation may be integrated once, and is then more conveniently expressed in 
terms of Q = R 2 : 

Q' 2 = 2(1 - fi')Q 3 + 4(/x'(Q) - 1)Q 2 + AHQ - 4ft 2 , (3.16) 

where H is a constant of integration. The solution is readily determined in terms 
of the Weierstrass elliptic function V (Abramowitz and Stegun 1965, Whittaker and 
Watson 1962) to be 

Q = -^T(X:g 2 ,g 3 ) + B, (3.17) 

where 

52 = -2(1- M ')^ + !(m'(Q>-1) 2 (3-18) 

.93 = 1(1 - m')(m'<Q) - m - 4f(M'(Q) - I) 3 + (1 - M') 2 ^ 2 (3-19) 

B= zWii), ( , 20) 

3(1 - M) 

and gi and 33 are the usual parameters of V . The amplitude and period of the solution 
are implicitly related through the appearance of the quantity (Q) in (3.18)-( |3.20 ). 



Although (3.17) provides an exact solution, it requires specification of //, H, ft 
and the spatial period to permit analysis of its solutions. Rather more analytical 
progress is possible if we restrict attention to the invariant subspace in which A is 
real. We suppose that the solution has spatial period L (note that, as a consequence, 



eL = 2n7r for some integer n). Equations ( J2 . 5| ) ( 2.6) become 

= A + Axx- A 3 ~AB, (3.21) 
= Bxx+v'(A 2 )xx, (3.22) 

where 

// = iij a. 



Integrating ( 3.22 ) twice and imposing periodicity gives the result 

B + fi'A 2 =f-i'{A 2 ), (3.23) 



corresponding to ( 3.14 ). Eliminating B from ( 3.21 ) gives 



= (1 - ti'{A 2 ))A + A xx - (1 - v')A 3 . (3.24) 
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Before describing the periodic solutions of (3.21)-( |3.22 ), it is instructive to 



consider a solution that is a good asymptotic approximation: 

A = A sech(cA), (3.25) 
where n' (A 2 ) = 1 + c 2 and fi' = 1 + 2c 2 jA\. Note that the solution exists only for 



[if > 1, which is consistent with the stability boundary ( pT9| ) and the result from § 3.2 
that the bifurcation at // = 1 is supercritical. For large values of cL, cL(A 2 ) ~ 2A^ 
Using this approximation for (A) , we obtain for c the quadratic equation 

c 2 - Afi'c/L{p' - 1) + 1 = 0, 

which has real solutions only for 

1<M' <L/(L-2), (3.26) 

indicating that a saddle-node bifurcation occurs at [i! = L/(L — 2) and that for 
/i' > L/(L — 2) no such solution exists. At large amplitude the solution branch 
approaches // = 1 and the amplitude is given by A) ~ 4y/2ii'/L{jj,' - 1) 3/2 . The 



solution (3.25) is of particular significance because, as we shall see, it can be stable 
and is observed in the numerical simulations described in § ^. 



Spatially periodic solutions of ( |3.24 ) can be written in terms of the Jacobi elliptic 



functions sn, cn and dn. In all that follows, to denotes the standard parameter of 
these functions (Whittaker and Watson 1962, Abramowitz and Stegun 1965), and we 
use the notation (• • •) as above, so that 



{f[X)) = ^f f(X)dX, 
where K is the real quarter-period, given by 

r' 2 de 



K(m) = 



(1-msin 2 0)V2" 



3.3.1. dn When fj/ > 1, there exist solutions of ( 3.24 ) in the form 

A = A dn cX, (3.27) 

where 

(2 - m)c 2 = Al(i'{ dncX) - 1, 2c 2 = A 2 Q (^ - 1), cL = 2K. (3.28) 



Given the parameter \j! and the length L of the domain, the three equations in (3.28) 
serve to determine the three unknowns A^ 1 c and to. However, the computation of 
solutions is more easily accomplished if L, say, is fixed, while to is varied in the range 
< to < 1 to yield corresponding values of A$, c and \i! . 
For small to, the solution takes the form 

A - (1 + \m){l - iTOsin 2 (7rA/ J L)), (3.29) 

indicating that these solutions branch from the uniform state A = 1; the branching 
takes place a t ji' = 1 + 2ir 2 /L 2 (corresponding to m = 0). Furthermore, it is readily 
shown from ( 3.28| ) that for small to 

M' = 1 + ^ + ^- 4 (L 2 - 4tt 2 )to 2 + 0(to 3 ), (3.30) 



confirming the result of § 3.2 that the bifurcation is supercritical when L > 2tt (I < 1 



in the notation of § [3.2;) and subcritical when L < 2tt (I > 1). 
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Figure 3. The bifurcation diagram for (3.21)— (3.22) with periodic boundary 
condit ions with L = 5ir (solid line) and L = IOtt (dashed line). Solutions take the 
form (3.27). The vertical axis shows the maximum value of A. 



As to is increased from zero to unity, the corresponding values of /z' increase to a 
maximum ^ (given approximately in (3.26) by [id ~ L/(L — 2)), and then decrease 
monotonically so that 

4 



1 



(3.31) 



ln(16/mi) 

as m — > 1~ , where mi = 1 — m. We have used the well known result that 

K ~ ln(4/m} /2 ) (3.32) 

in this limit. Thus a sad dle-n ode bifurcation takes place at fi' — Hd, and for // > /id 
no solutions of th e form ( 3.27 ) exist. Figure || shows the bifurcation diagram for dn- 
type solutions of ( 3.21 )-( p^22 ) subject to periodic boundary conditions with L = 
and L = IOtt. The saddle-node bifurcations occur at /id = 1.146 and fid = 1-068 
respectively, in agreement with our approximation (3.26). Since the bifurcation is 
supercritical in each case, the lower branch of solutions is stable but the upper branch 
is unstable. 

When m is small, the solution for A takes the form of small, almost sinusoidal 
perturbations to the uniform state A = 1. As to is increased, the peak of A becomes 



narrow and the trough broadens. As to 

and the solution takes the form 

„ 2K 3 ? 2 , 2KX 
A ~ — - — sech — = — 



1 , the peak increases greatly in amplitude 



(3.33) 

In this limit, the solution for 



where [i! ~ 1 + 2/K, and K — > oo according to ( [3.3 
A thus takes the form of a single sech-type pulse as given by (3.25). The dn solutions 
are illustrated in figure |[ 



These results indicate that for > // > 1 a soliton-like solution exists to (2.5)- 
( [2.6| ). For /i' > /id, numerical simulations of (2.5)— (2.6) indicate that the solution 
blows up, signalling a transition to a fully nonlinear solution of the original PDEs 
(for example (2.7)). This blow-up can be avoided by adding higher-order terms to the 
amplitude equations, which we discuss in § 0. 
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Figure 4. Solution of (3.21) and (3.22) in terms of dn (3.27) for m = 0.1, 0.5, 0.9, 
0.999, 0.9999999 (corresponding to~if= 1.0200, 1.0203, 1.0230, 1.0504, 1.0596). 
The spatial period is L = 31.4. 




Figure 5. The function g(m), defined in 
given by g(0) = 4tt 2 and g(mo) = 0, where mo 



The intercepts with the axes are 
s 0.826115. 



3.3.2. cn When \J > 1, solutions also exist in the form 
A = A cn cX, 

where 



(2m - l)c 2 = A%(i'{ cn 2 cX) - 1, 
When m is small, A ~ v2(l 



2mc 2 = Al(n' - 1), 



cL = AK. 



(3.34) 



(3.35) 



47T 2 /!/ 2 ) 1 / 2 cos(27rX/L), and hence we require 
L > 27r for the existence of this solution in this limit. Further analysis indicates that 
for general values of m a solution exists only when 

L? > g(m) = \QK 2 [l + 2m(( cn 2 cX) - 1)]. (3.36) 

The function g{m) is shown in figure |5|. It is clear from ( [3.36 ) and figure || that 
if L > 2-7T the solution ( p. 34 ) exists for all m between zero and unity. If, on the 
other hand, L < 2tt then solutions exist for mj, < m < 1, where g(rriL) = L? . 
Regardless of the value of L, these solutions exist for m > mo ss 0.826115. If £ > 2ir 
then as m is increased from zero, // increases from unity to a maximum value [i c 
(given approximately by fi c = L/(L — 4)) and then decreases monotonically according 




Figure 6. Solution of (3.21) and (3.22) in terms of cn (| 
0.999, 0.999999 (corresponding to /7 
period is L = 31.4. 



^34|) for m = 0.1, 0.9, 
I.. Out I.. 0(i9. 1.1.28,1.146). The spatial 



to (3.31) as m — ► 1 . If L < 2tt then as to is increased from m^, // decreases 



monotonically from infinity, again according to (3.31) as m — > 1~. 

As to is increased, the peaks of A become more pronounced, larger in amplitude 
and narrower, while the regions in which A has small amplitude become broader. In 
the limit m-tl", the solution is approximately given by two sech-type pulses, one 
positive and one negative; with an appropriate choice of origin in X, the solution for 
< X < L satisfies 

4^3/2 



A 



, AKX , AK(X-L/2) 
sech — sech 



with K and m related through ( 3.32| ). The cn solutions are illustrated in figure |(| 



3.3.3. sn The final type of solution that may be expressed in the form of Jacobi 
elliptic functions applies when fx' < 1, and is 

A = A sncX, (3.37) 

where 

(l + m)c 2 = l-A&'i sn 2 cX), 2mc 2 = A 2 {1 - fj,'), cL = AK. (3.38) 

For small m, these solutions match with the cn solutions (allowing for an 
unimportant translation in X), with A ~ v2(l — 47r 2 /L 2 ) 1 / 2 sin(27rX/L). As this 
expression indicates, the small-m solutions require L > 2ir. In fact, for general values 
of to, the length L of the domain must satisfy 

L 2 > h(m) = 16K 2 [1 + to - 2m( sn 2 cX)}. (3.39) 

The function h(m) is monotonic in creas ing, with h(Q) = Air 2 and h(m) — > oo as 



to — > 1 . Thus solutions of the form ( 3.37 ) do not exist if L < 2tt; if L > 2n, solutions 



exist for < to < tol, where /i(tol) = L . As to increases towards tol, fi' decreases 
monotonically to — oo, the peaks and troughs in A becoming broader. Provided L is 
large, the profile near X — is approximately a front of the form 

, 4KX 
A ~ tanh ■ 



L 

the front near X — L/2 being given by symmetry. The sn solutions are 
illustrated in figure [j]. 
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Figure 8. The bifurcation diagram for (3.21 )—( 3.22 ) with periodic boundary 
conditions with L = IOtt, showing dn, cn (both solid lines) and sn (dashed line) 
solution branches. The vertical axis shows the maximum value of A. 



A bifurcation diagram showing all three types of solution (cn, dn and sn) is given 
in figure |^. Note that only the dn solutions bifurcate from the periodic rolls represented 
byA=l. 

The solutions described above are those containing the minimal number of peaks 
and troughs within the periodic domain. Other solutions, containing n sets of peaks 
and troughs can also be constructed by replacing K by nK in the last equation of 
( p8| ), fU5j ) or d3~38) ). 



4. Numerical simulations of a model PDE 



In this section we illustrate the application of the results of the previous sections 
through numerical simulations of the PDE (2/7), assuming periodic boundary 
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Figure 9. Numerical simulations of (2.7) showing w(x), with r = 0.01 and 
periodic domain size 62tt. (a) Stable rolls, s = 0.9; (b) Modulated rolls, s = 0.93; 
(c) s = 0.95; (d) s = 1.0. 



conditions. From the analysis of § 2.2, the system is governed by the amplitude 



equations (2.5)— (2.6) with 
1 



and 



^ = 



(4.1) 



3 - 2s 2 /9 

A pseudospectral method is employed in which the nonlinear terms are computed in 
physical space using FFTs. Time advancement is carried out by approximating the 
nonlinear terms as a linear function of time and then obtaining exact solutions using 
the known growth rate of the linear terms (2.8). 

In order to approach the asymptotic regime, the size of the periodic domain must 
be large and the driving parameter r must be small. Initially, the length of the domain 
is chosen to be 627r, so that 31 periods of the pattern are contained in the box, the 
forcing parameter is r — 0.01 and the number of grid points in physi cal s pace is 256. 
The effective length of the domain for the amplitude equations ([2.5|)-(2.6) in terms of 
the leng thscale X is then L = 62iry/r/2 — 3.1tt. Note that since L > 2ir, the result 
of § 3.2 indicates that the bifurcation from a regular pattern to a modulated pattern 
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is supercritical. In the limit of large L the condition (3.9) for rolls to be unstable is 
/j, > cr, which is equivalent to s 2 > 27/38. For finite L the condition is (3.6), which 
yields /i > er(l + 2ir 2 /L 2 ), or |s| > 0.9215 for L — 3.1ir. In our numerical simulations 
of (2.7) at r = 0.01, we find good agreement with this asymptotic result, rolls being 
unstable for |s| > 0.927. 

A sequence of numerical solutions of (2.7) with r — 0.01 for different values of 
s is shown in figure ^. In all cases, the solution shown is the final steady state after 
all transients have decayed; to reach this state the computation must be continued 
for several thousand time units. In figure |§(a), s — 0.9 and a stable periodic solution 
is found. For s = 0.93 (figure ||(b)), slightly beyond the bifurcation, a modulated 
pattern appears. The amplitude of the modulation is small, which is consistent with 



the result of § 3.2 that the bifurcation is supercritical. As s is increased, the amplitude 



of the modulation steadily increases (figure |9j(c,d)), until there are regions in which 
the pattern is almost completely suppressed and regions in which the amplitude of 
the pattern becomes large. The envelope of the pattern takes the form of the dn or 
sech solutions described in § |3.3| (see figure ||) ; we have not found any of the cn or sn 
solutions, which suggests that they are unstable. 

The increasing amplitude of these solutions means that the numerical simulations 
arc no longer within the asymptotic regime of validity of the equations ( |2.5[ )-(2.6); 
in particular, for this value of r we do not find the saddle-node bifurcation, which, 
according to (3.26), should occur when s — 0.939 in (2.7). 

To find this saddle-node bifurcation it is necessary to carry out further sequences 
of numerical simulations at smaller values of r. In order to keep the value of L 
constant, the length of the domain (and the number of grid points) must be increased 
in proportion to r -1 / 2 . The numerically obtained bifurcation diagrams, with r 
successively reduced by a factor of 4, are compared with the solution of the amplitude 
equations in figure [l^. The function A is obtained numerically from the part of the 
spectrum of the numerical solution near k = 1 and scaled appropriately with r 1 / 2 . For 
a domain size 62-7T with r = 0.01 there is no saddle-node bifurcation and the numerical 
and theoretical bifurcation diagrams are not in close agreement. With a domain size 
1247T and r — 0.0025, agreement is better near the pitchfork bifurcation but there is 
still no saddle-node bifurcation. Finally, for a domain size 2487T with r = 0.000625 a 
saddle-node bifurcation appears at s — 0.944, close to the theoretical value of 0.939. 
Beyond this point the numerical results jump to a solution of much larger amplitude, 
and there is a region of hysteresis where two stable sol ution s exist. 

For larger values of s, finite- amplitude solutions of (2.7) are found, that cannot be 
described in terms of the amplitude equations. In some cases there are multiple stable 
solutions. For example, in figure [ll] we show two solutions computed at s = 1.6, from 
different initial states. In one plot there are two regions of pattern, and in the other 
only one. The amplitude of the solution appears to be approximately the same in 
each simulation. However, the single pulse has roughly twice the width of each of the 
smaller pulses, so in each simulation an equal proportion of the computation domain 
seems to be occupied by the pattern. Furthermore, it is clear that the envelope is 
rather broader than a sech profile, particularly in figure |ll|(b), where the envelope 
consists of two nearly uniform regions separated by a pair of fronts. In addition to 
those displayed in figure [n], we have found further steady (or quasi-steady) solutions 
of (2.7), including a single-pulse solution containing six large waves (rather than the 
seven in figure |ll](b)) and a two-pulse solution where the pulses are of unequal size, one 
containing three, and the other four large waves. Closer examination of the numerical 
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(a) 

1 



(b) 




-2 



Figure 11. Numerical simulations of (|2J) showing w(x), with r = 0.01, domain 
size 62-7T (using 512 grid points) and s = 1.6. The two solutions are started from 
different initial conditions. The single-pulse solution appears to be a steady state, 
while in (a) the two pulses are slowly separating. The amplitude appears to be 
almost uniform in each of the humps. 
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results indicates that the multi-pulse solutions are not quite stationary; the pulses are 
gradually moving. 

5. Beyond the leading-order amplitude equations 

In the previous section we have seen that agreement between the amplitude equations 



( |2 . 5[ ) — ( J2 . 6| ) and the behaviour of the full PDE (2.7) only occurs for extremely small 
valu es of r . This agreement can be improved by introducing higher-order terms into 
( |2.5|) (|2.6|). R ather than include all possible terms that might arise at the next order 
in ( |2.5| )-( |2.6[) , we focus on one term. 

To motivate our choice of higher-order term, we note that when we compute the 
x- average of w 3 , the leading-order contribution is e 4 (|s|^4| 4 + 6|A| 2 i?). For moderate 
values of s, the term 6e 4 |A| 2 i? dominates. After rescaling the equations for A and B 



so that A satisfies (2.5), we find that B satisfies 

B T = aB xx + ti{\A\ 2 ) xx + 8{B\A\ 2 ) XX , (5.1) 



where er and fi are given by (4.1) and 



3f 2 

5=—. ^v>0. (5.2) 

2(3- § s 2 ) 

For future convenience we introduce the quantity M, defined by 

((M-^'\A\ 2 )/(l + 6\A\ 2 /a))=0. (5.3) 

When \j! > 0, the term in 5 plays a stabilising role. 

For steady solutions, the governing equation for A becomes 

O-A + Asz-AW-A^M^W). (5.4) 



The bifurcation diagram for this equation, subject to the constraint (5.3), is indicated 
in figure [l2| The bifurcation diagram is changed significantly even by very small values 
of 5. For example, whereas there is a single saddle-node bifurcation for 6 = 0, when 
S > there is a pair of saddle-node bifurcations. This pair exists only for 5/a less 
than approximately 0.0009, and for larger values of S/a the curve is monotonic. 

We now wish to determine the effect of small 5 upon the envelope of the pattern, 
since our numerical simulations suggest that the envelope will broaden from a sech 
profile as the envelope becomes more developed. For small 5, J5.4| ) is approximately 

= (1 - M)A + A xx - (1 - // - SM/a)A\A\ 2 - ^'8A\A\ A /a. (5.5) 

and so to determine the effect of 5 upon the envelope, it is instructive to examine the 
model equation 

f" = f-f + Sf. (5.6) 

This model contains the essential features of the problem ( |5.3| ) and (5.5) for A, and 
has the exact solution 

f = b 1 / 2 (}(b-2) + cosh 2 xy 1/2 , (5.7) 

where b = 2/(1 — 16J/3) 1 / 2 . When 5 = 0, / = \/2 sechx, but as 5 increases the 
profile broadens and flattens (see figure Q), giving envelopes reminiscent of those in 
figure [n]. 
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Figure 12. Bifurcation diagram for (Q subject to for 5/a = 10" 4 (solid 

line), 10~ 3 (dashed line) and 10" 2 (dotted line) plotted using AUT097. The 
domain length is L = 5n. Note that very small values of <5 can change the 
qualitative nature of the bifurcation diagram, for example eliminating the saddle- 
node bifurcation at 8 = 0. 




Figure 13. Envelopes ( p.7[ ) from the model equation (5.6) with 5 = (solid line) 
and S = 0.185 (dashed line). The larger value of <5 produces a broad profile with 
a nearly uniform region inwhich / is non-zero. The broader profile is reminiscent 
of the envelope in figure hlL 



6. Examples 

Since conservation laws arise very commonly in nature, it turns out that a wide variety 
of pattern formation problems are governed by the system (2.5)— ( pT^ ) rather than the 
Ginzburg-Landau equation alone. 

In thermal convection in a rotating layer with stress-free boundaries, the 
conserved quantity is horizontal momentum and the corresponding long-wave mode 
is a large-scale flow along the axes of the convection rolls. In the analysis of Cox 
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and Matthews (2000), ( |2.5| )-( |2.6| ) appear as a special case, and it can be shown that 
rotating convection is unstable for moderate values of the Prandtl number. 

In Marangoni convection (driven by surface tension), the fluid has a free 
surface and the equation for the layer depth is constrained by conservation of mass. 
Equations similar to ( |2.5| )-( f2~6| ), but with one extra term, were derived by Golovin, 
Nepomnashchy and Pismen (1994); in this case, however, the parameter (i is negative 
and so no instability arises. 

In §§ 6.1 an d |6.2| below, we derive an equation similar to ( |2.7| ) for two p robl e ms in 
convection; in § p.3| , we directly derive amplitude equations equivalent to (2.5)— (2.6) 
for magnetoconvection between 'ideal' horizontal boundaries. 



6.1. Application to convection with a nonlinear equation of state 

A physical problem that gives rise to amplitude equations of the form considered 
here is convection between boundaries at which the heat flux is fixed (Roberts 1985, 
Matthews 1988a, 1988b). The fluid has a nonlinear equation of state, 

p. = Po [1 - a(T, - T ) - 0(T, - T ) 2 ] , (6.1) 

relating the fluid density to the temperature T*. The coefficients a, f3 and T are 
constants. This equation of state is appropriate for water around 4°C, for example. 
The crucial feature of this problem is that the total heat in the layer is conserved; 
the fixed-flux boundary conditions allow a large-scale, slowly-evolving depth-averaged 
temperature mode, 0. 

We suppose that the temperature profile is nonlinear, of the form (Matthews 
1988a, 1988b): 

T* = T + Az* - Bz*, (6.2) 

where z* is the vertical coordinate. The problem is expressed in dimensionless variables 
by adopting the scales d = (B/A) 1 ^ 2 , Bd and d 2 /n for length, temperature and time. 
We then find the equations governing two-dimensional motions to be 

P- 1 (VV* + ^(V-,V 2 ^)) = i?^ + V> + i?A(z 3 -z + 6)6 x , (6.3) 

t + J(<AA) + (3z 2 -l)^ = V 2 #, (6.4) 

where the x- and z-components of velocity may be expressed in terms of the 
streamfunction ip as (it, w) — (—ip z ,ip x ), and the Jacobian is defined by J(f,g) = 
fx9z — fz9x- The perturbation to the temperature profile ( |6.2[ ) is denoted by 9. The 
Rayleigh number is 

gaBd 4 
H = , 

where v and k are, respectively, the coefficient of kinematic viscosity and the coefficient 
of thermal diffusivity, g is the gravitational acceleration; the Prandtl number is 
P = v/k, and A = 2Bd(3/a. In deriving the PDE for Q, we assume that the equation 
of state is only mildly nonlinear, so that | A| <C 1. 

We break the up-down symmetry of the problem by imposing rigid boundary 
conditions at z — —H but stress-free conditions at z = H. Thus 

j/j = i/} z = at z = -H (6.5) 

and 



ip = ijj zz = at z = H. 



(6.6) 




Figure 14. Marginal stability curves for convection with a nonlinear equation 
of state and nonlinear basic-state temperature profile. When the depth H of the 
fluid layer is less than some threshold Ho, the critical wavenumber is zero and 
the critical Rayleigh number is Ro, given by ( p,7| ); when H > Ho, the critical 
wavenumber k± is positive, and the critical Rayleigh number i?i is less than Ro- 



The important mathematical consequence of these boundary conditions is that the 
equation for then contains a quadratic nonlinear term of the form {Q 2 x )xx, where 
X is a scaled ir-coordinate. This term tends to lead to a supercritical bifurcation to 
rolls, unlike the quadratic term (9 2 )xx, which arises from the nonlinear equation of 
state, and which tends to lead to a subcritical bifurcation. 

We suppose that the heat flux through the planes z — ±H is fixed, so that the 
perturbation heat flux 89 /dz — at each boundary. The key consequence of this 
choice of boundary condition is that, for sufficiently small half-depth H, the critical 
mode of the linear stability problem has zero wavenumber (see figure |lj) . 

We begin by considering the si mpl er case of a linear equation of state, with A = 0. 
First we examine equations (5.3)-( [6~4| ), linearised about the trivial state ip = 6 = 0. 
We find, making asymptotic expansions of the forms given by Roberts (1985) and 
Matthews (1988a), that modes with wavenumber k have linear growth rate A, where 

A ~ (R - R^R^k 2 - ai k 4 , Ro = 420iJ~ 4 (lliJ 2 - 21) _1 , (6.7) 
for small \R — Rq\ and k. The coefficient 04 is 

_ 8£T 2 (6831iJ 4 - 25738£T 2 + 23751) 
° 4 ~~ 1287(lli? 2 - 21) 2 ' 

As a consistency check on these results, we note that as H — ► 0, Ro ~ 20/ H 4 and 
ct4 ~ 232i7 2 /693. This limit corresponds to the problem with a linear equation of state 
and linear basic-state temperature profile, such as considered by Chapman and Proctor 
(1980). After appropriately rescaling of the depth of our fluid layer, we find that their 
results agree with ours in the limit as H —> 0. When 04 > 0, the largest wavenumbers 
in (3.7) are damped, and k = corresponds to a minimum of the marginal stability 
curve. This condition is satisfied for < H < H , where H w 1.2709 is the smaller 
positive root of 04 = (see figure 14). 

The amplitude equation of interest applies when H is close to H . We thus follow 
Roberts (1985) and Matthews (1988a) and let 

R = R + e 2 R 2 + e 4 Ri, A = e 2 D, H = H a + e 2 H 2 , (6.8) 
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expanding the streamfunction and temperature perturbation as 
8 =e 2 [8 2 (X, T) + e 2 6 4 {X, z, T) + ■ • •] 

yj 3 (X,z,T) + e 2 MX,z,T) + ---] , 



(6.9) 
(6.10) 



where X = ex and T — e 6 t. These scalings have been adopted to ensure a balance 
of terms in the equation derived below. Of course, it is reasonable to question the 
physical relevance of such a large power of e in the scaling of time, and we offer two 
partial justifications. The first is that, in fact, solutions to model PDEs of the form 
to be derived below are known to provide a good approximation to the full convection 
problem even when e is not particularly small (and hence e 6 represents a physically 
achievabl e sca ling). The second justification is that in deriving amplitude equations of 
the form (2.5)-( [2.q ) it is not necessary to pass through an intermediate stage involving 
the time scale e 6 t, this is just a convenient way to construct in an asymptotic fashion 
a long- wavelength mo del PDE. 

At 0{e 2 ) in |y|-(fl]), we find that d 2 2 /dz 2 = 0, subject to the boundary 
conditions 882 /dz = at z = ±H. The solution for 82 is therefore independent of 
depth, as indicated in (|6.9D . After some algebra, we find from a solvability condition 
at 0(e 8 ) that the amplitude equation for = 9 2 (X,T) is 

^ + a 2 D 2 ) 6 + [a' A H 2 + a'{D)Q xx + a 6 Q X xxx + a q G 2 x - ±DQ 2 \ 
u o J ) XX 

where the numerical values of the coefficients in ( |S.ll ) are given approximately by 

a 2 = 7.588 x 10~ 5 , (6.12) 

a' A = 8.966, (6.13) 

04 = 0.1208, (6.14) 

a 6 = 0.3067, (6.15) 

a q = 5.571 - 1.086/P, (6.16) 

Rq = 49.80. (6.17) 

To ensure that all terms in the equation for Qt arise at the same asymptotic order, 
i?2 = CL7D, where 07 w 0.4338. 



e 



(6 



After rescaling, we can write (6.11) as 



dw 
~dt 



dx 2 



1 



dx 2 



w — sw — 



dw 
dx 



(6.18) 



which is very similar in form to (2.7) considered in § 2.2. Following the asymptotic 
procedure outlined in § 2.1 , we find that with w(x,t) as in ( [2.10 ), the amplitudes of 
the pattern and large-scale modes satisfy 

A T = r 2 A + AA X x - 2(1 - s)(s + 2)A\A\ 2 /9 - 2sAB, (6.19) 

(6.20) 



B T = B xx + 2{s + l){\A\ 2 ) xx . 



If —2 < s < 1, the coefficient of the term A|yl| 2 in (6.19) is negative, and the 
bifurcation to pattern (in the absence of B) is supercritical. Equations ( |6.1£ ) and 
( |6.20 ) are then equivalent to (2.5)— (2.6), after an appropriate rescaling. 

From this analysis we find that the bifurcation to rolls is supercritical if — 2 < 
s < 1, and that all rolls are unstable if 



2 < s < 



3V57 
38 



-1.096 or 0.096 : 



3V57 
38 



< s < 1, 



(6.21) 
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where 

1/2 R 1/2 A 

2a q (R + a 7 A + a 2 R A2~R) 1 / 2 ' 1 ' ' 

6.2. Application to thermosolutal convection 

As a second physical example, we examine thermosolutal convection between 
boundaries at which the heat flux is fixed. We suppose that the fluid is slightly non- 
Boussinesq, so that the thermal diffusivity is a weak function of the local temperature: 

k* = k (1 + k(T*-T )/T ). (6.23) 

The layer is heated from below and in the basic state the saltier water lies over the 
fresher water, so both mechanisms are destabilising. Both the thermal and the solutal 
gradients are vertical and uniform in the absence of convection. The thermal gradient 
is characterised by a (thermal) Rayleigh number R > and the salinity gradient by 
a solutal Rayleigh number S < 0. When > S > S c , the bifurcation to convection 
takes place at zero wavenumber; when S c > S, the critical wavenumber is finite. We 
analyse the case 

S = S c + e 2 S 2 . 

The natural scale for lengths is now the depth of the layer, and that for 
temperature perturbations is the temperature difference across the layer in the basic 
state (similarly for the salt concentration). With these scales, the dimensionless 
governing equations may be written in the form 

P- 1 (vVt + J(iJ>, V 2 ^)) = R9 X - SV X + VV (6-24) 

0t + J(V>, 0) - V* = V 2 6» + kV • [(-z + 0)V0] (6.25) 

E t + J(V>, E) - V* =rV 2 E, (6.26) 

where E is the perturbation salinity field, r = k s /k and k s is the salt diffusivity. The 
thermal Rayleigh number is now 

R= ga(dT*/dz*)d^ 

where a is the coefficient of cubical expansion and dT^/dz* is the (dimensional) 
temperature gradient in the absence of convection; the solutal Rayleigh number S 
is defined similarly. 

As in the previous example, we use free-rigid horizontal boundary conditions, with 
( |6.5[ ) and ( |6.6| ) applied for H = 1/2, since the fluid layer now has unit dimensionless 
depth. With these boundary conditions, S c = -8352r/41 m -203.707r. 

The amplitude equation of interest applies when 

R = i?o + £ A Ra, S = S c + e 2 S2, k = e 2 K,2] 
we expand the independent variables as 

9 = e 2 [0 2 (X, T) + e 2 9 4 (X 7 z,T)-\ ] (6.28) 

E = e 4 [E 4 (A, z, T) + e 2 E 6 (A, z, T) + ■ ■ •] (6.29) 
rl> = e 3 [MX, z, T) + e 2 ^ 5 (A, z,T) + ■ ■ ■] , (6.30) 
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where X = ex and T = e 6 t. At 0(e 8 ) in an asymptotic expansion of ( 6.24 )— ( |6.26 ) , we 
find that the amplitude equation for = 9 2 (X, T) is 



Gn 



Ro 



+ c 4 S 2 &xx + c e Qxxxx + c q 2 x + \k 2 & 



.(6.31) 



XX 



where 



Ro = 320, 



41 



c 4 



-0.0004109/r, 
0.01312, 



99792r 
14506722953 

1105865601750 
5 4- 7P 

« -0.08333 - 0.05952/F. 

84P ' 



(6.32) 
(6.33) 

(6.34) 

(6.35) 



Note that this has exactly the same form as (6.18). From § 6.1, the condition for a 
supercritical bifurcation is — 2 < s < 1, and that for all rolls to be unstable is (6.21), 
where 



2c q (R -R) 



i/2- 



(6.36) 



6.3. Application to magnetoconvection 

Convection in a vertical magnetic field has been a subject of a great deal of research 
(see, for example, Proctor and Weiss 1982; Matthews 1999). The problem is motivated 
by the convection zone of the Sun, where the magnetic field is moved around by the 
fluid and regions of strong magnetic field inhibit the fluid motion. In the nonlinear 
regime magnetoconvection exhibits a very wide range of complicated phenomena. 

In magnetoconvection, the conservation law that leads to the large-scale mode is 
that the total flux of magnetic field through the layer is conserved. The governing 
equations for two-dimensional magnetoconvection are (Proctor and Weiss 1982), in 
the notation of the preceding section, 

P- 1 (V 2 ^ + J(^VV)) =i?^ + vV + CQ(^V 2 0)+V 2 z ), (6.37) 

6 t + J(iP,6) =iP x + V 2 6, (6.38) 

<j)t + J(ip, <t>) = i>z + CV 2 ^. (6.39) 

The magnetic field is defined through a perturbation flux function 4> by (B X ,B Z ) = 
(— 4> z ,l + 4> x ). The parameter Q is the Chandrasekhar number, a dimensionless 
measure of the magnetic field strength, while £ is the ratio of magnetic diffusivity 
to thermal diffusivity. 

The most commonly used boundary conditions are that the upper and lower 
boundaries are stress free, fixed temperature and with no horizontal magnetic field, so 

i\> = i\) zz = 6 = <j) z = at z = 0, 1. (6.40) 

These boundary conditions are convenient in the sense that the linear eigenfunctions 
are trigonometric. 

The large-scale mode in the problem arises through a rearrangement of the vertical 
magnetic field. There is a linear eigenfunction in which <f> — <f>(x, t) which obeys the 
diffusion equation (f>t — C4>xx , so the growth rate A of a mode with wavenumber k is 
simply A = — (k 2 . Note that this mode is independent of the velocity and temperature 
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perturbation, since if <fi — fax,t) there is no forcing term in (6.37). Although this 
mode is usually ignored in the analysis of magnetoconvection, it plays a crucial role 
when the horizontal size of the domain is large. Note that although there are nearly 
marginal modes at small fc and at order-one values of fc, the graph of growth rate 
versus wavenumber is not as shown in figure |l], because the large-scale mode and the 
order-one mode correspond to two distinct branches of eigenfunctions. 

The weakly nonlinear analysis proceeds by including the large-scale mode. The 
general analysis of §2 suggests that the mean mode should be included at order e 2 . 
In fact the mean vertical magnetic field is of order e 2 and so the corresponding flux 
function cj) is of order e. The appropriate asymptotic expansions near the stationary 
bifurcation to convection are therefore 

tp = eA(X,T) exp ikx sin ttz + c.c. + 0(e 2 ), 

9 = eA(X,T) 6>iiexpifcrshi7rz + c.c. + 0(e 2 ), 

cj) = eA(X, T) fa exp ikx cos nz + c.c. + e fa T) + 0(e 2 ) 

where 6i and fa are constants, readily determined from the linear stability problem, 
X = ex, T = e 2 t and R = R c + e 2 R 2 . At order e we obtain 6 X = k/(n 2 + k 2 ), 
(pi = tt/((tt 2 + fc 2 ) and the well known result (Chandrasekhar 1961) 

(n 2 + k 2 ) 3 + Qir 2 (iT 2 + k 2 ) 



(6.41) 
(6.42) 
(6.43) 



Rc 



(6.44) 



From ( |6.44 ) it follows that as R is increased, convection first occurs with a critical 
wavenumber k c which is related to Q by 

Qtt 4 = (2fc 2 - 7r 2 )(7r 2 + fc 2 ) 2 . (6.45) 

Henceforth it is assumed that k = k c ; note that the envelope function A(X, T) permits 
convection with a wavenumber k — k c + 0(e). 

At order e 2 , a term in proportional to sin27rz and a <fi term proportional to 
exp 2ikx are obtained. At third order, evolution equations for A and $ are found by 
applying solvability conditions. Full details of the calculation will be given elsewhere. 
The resulting amplitude equations are 

A T = a x A + a 2 A X x - a 3 A\A\ 2 - a 4 A$ x , (6.46) 

$t =C<5>xx+k(\A\ 2 ) x , (6.47) 

where a% . . . 0,4 are known constants, given by cumbersome functions of the parameters; 
02 and a 4 are always positive. After setting $x = B, ( |6.46 )-(6.47) are equivalent 
to (2.5)-(fOi|). The equations (6.46)— ( |6.47 ) have not been derived previously, but the 
underlying processes of the coupling between the convection rolls A and the large-scale 
mode $ are well unde rstoo d (see, for example, Proctor and Weiss 1982). Since 04 > 0, 
the coupling term in (6.46) simply reflects the fact that the magnetic field suppresses 
convection. The coupling term in (6.47) describes the process of flux expulsion: the 
magnetic field strength decreases in regions where the convection is strong. 

The bifurcation at the onset of convection is supercritical (i.e. 03 > 0) provided 

that 

fc 4 CV + k 2 ) > 7T 2 (2fc 2 - 7T 2 )(7T 2 - fc 2 ), (6.48) 



where Q has been eliminated in favour of k using (6.45). The result (3.9) can now be 
applied to determine whether a regular pattern of rolls is unstable. It is found that 
all rolls are unstable if 

fc 4 CV + fc 2 ) < 7r 2 (2fc 2 - 7r 2 )(fc 2 + 3^ 2 ), (6.49) 
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Figure 15. Region of stability of convection rolls in magnetoconvection. To the 
right of the solid line, rolls are stable. To the left of the dashed line the bifurcation 
is subcritical. Between the two lines, all rolls are unstable. 



where Q has again been eliminated using ( 6.45 ). Note that neither of the conditions 
( |0.48 ) -( 6.49| ) involves the Prandtl number P. Using k as a parameter, the region 
of stability of rolls can be plotted in the Q, ( plane (figure 15). In the small region 
where the rolls are subcritical the weakly nonlinear analysis gives no useful information 
about the pattern. However, there is a large region of parameter space in which the 
bifurcation is supercritical but all rolls are unstable. In this region we expect to find 
that the stable pattern at onset consists of rolls that are modulated on a large length- 
scale. Solutions of this type, but fully nonlinear, have recently been obtained through 
numerical simulations of magnetoconvection (Blanchflower 1999). 

Note that our analysis is confined to the onset of convection at a stationary 
bifurcation. For certain parameter values (in general, when Q is large and £ is small) 
the onset of magnetoconvection is oscillatory; since this depends on the value of P, 
this condition is not shown in figure [l^. 



7. Discussion 



We have derived amplitude equations (2 : 5)-(^_^) for the primary bifurcation to a 
stationary pattern in a one-dimensional system with a conservation law for a quantity 
that does not change sign under reflection. In this case the conservation law leads to a 
neutral large-scale mode which must be taken into account for a correct description of 
the behaviour near onset. The usual Ginzburg-Landau equation (Melbourne 1999) is 
coupled to an equation describing the large-scale evolution of the conserved quantity. 
We have shown how these equations may be obtained from symmetry considerations 
alone, and have also derived them directly from the governing equation(s) of the 
system. 

In considering the small-amplitude solutions of the amplitude equations, we have 
shown how periodic patterns may become unstable to a new type of instability. This 
new instability is quite different from the usual Eckhaus instability in three ways. It 
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is an amplitude-driven instability (whereas the Eckhaus instability is phase-driven) 
and it is in general supercritical. Most importantly, for certain parameter values 
all patterns are unstable. Our instability is also different from that which occurs in 
the case where the conserved quantity is velocity-like, and the reflection symmetry 
acts as —1 (such as in systems with Galilean invariance), where periodic patterns are 
always unstable and chaotic dynamics is observed at onset (Tribelsky and Velarde 
1996, Tribelsky and Tsuboi 1996, Matthews and Cox 2000). 

Beyond the instability, there are stable solutions in the form of a strongly 
amplitude-modulated periodic pattern, the envelope of which we have calculated in 
terms of Jacobi elliptic functions. The envelope may commonly be approximated by a 
sech profile. Numerical simulations reveal that these strongly modulated patterns 
may be stable in the full governing PDE, not just within the framework of the 
amplitude equations. Numerical simulations of the governing PDE also demonstrate 
that the amplitude equations represent a good approximation provided the amplitude 
is sufficiently small, but that higher-order terms play an important role as the 
amplitude increases; these higher-order terms can lead to significant modifications 
in the bifurcation diagram, and a broadening of the envelope of the pattern. 

Our work demonstrates that localised solutions are possible when a Ginzburg- 
Landau equation is coupled to an equation for a mean field, even when the coefficients 
of the equations are real and when the bifurcation is supercritical. Previous work has 
required complex coefficients (Riecke 1992a, 1992b, Sakaguchi 1993), a subcritical 
bifurcation (Herrero and Riecke 1995, Sakaguchi and Brand 1996) or both. Our 
amplitude equations apply when the bifurcation from the basic state is stationary; 
the case of oscillatory bifurcation to travelling-wave solutions has been analysed in a 
series of papers by Riecke and co-workers (Riecke 1992a, 1992b, 1996, Herrero and 
Riecke 1995) and by others (Sakaguchi 1993, Barthelet and Charru 1998), who find a 
variety of localised and front-type solutions. They find that the coupling to a mean 
field (our B) can have a significant effect on solutions and their stability, and can 
stabilise localised travelling-wave trains, and holes. 

We have applied our amplitude equations to three convection problems, leading to 
novel instabilities in these systems. Further details of these applications will be given 
elsewhere. In the first two, the heat flux is fixed on the boundaries, and the total heat 
in the layer is conserved, which leads to a large-scale mode in the form of a depth- 
averaged temperature perturbation. However, in order for our amplitude equations 
(|2.5|)- (|2.6| ) to apply, this large-scale mode must influence the convection rolls; this 
is achieved by using either a nonlinear equation of state or a temperature-dependent 
diffusivity. In the third example, magnetoconvection, the conserved quantity is the 
vertical magnetic field; in this case the large-scale mode and the pattern mode 
correspond to two distinct branches of eigenfunctions of the linear problem. Our 
modulated solutions may be related to the isolated 'convectons' found recently in 
magnetoconvection (Blanchflower 1999). Another application of our work is in rotating 
convection, where the conserved quantity is momentum (Cox and Matthews 2000). 

Our amplitude equations fl2.5|)-([2~6|) are identical to those derived by Coullet 
and Iooss (1990), who considered the secondary stability of a one-dimensional cellular 
pattern. In that case, A corresponds to the amplitude of a perturbation to the basic 
pattern, and B to a gradient in the phase of the pattern; the amplitude equations 
then describe evolution near a stationary bifurcation. Coullet and Iooss (1990) 
found that when a pattern of wavelength L becomes unstable to a perturbation of 
wavelength 2L, the amplitude equations (2.5)-([2~6|) apply, with A real; when the 
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spatial period of the perturbation and the basic pattern are irrationally related, A 
may be complex. Our results suggest that if a pattern of wavelength L undergoes 
a supercritical subharmonic bifurcation then there may be stable states where the 
original pattern is largely undisturbed in some regions but where the perturbation is 
strong in others. Furthermore, provided the secondary bifurcation is supercritical, it 
should be possible to find stable states consisting of two domain types: in one domain 
the pattern is essentially the basic cellular state, in the other the pattern is strongly 
perturbed. 

It would be of interest to investigate how this theory may be extended to problems 
in two or three space dimensions; in particular we would like to know whether solutions 
may be obtained that are localised in space. Such localised solutions have been found 
experimentally by VanHook et al. (1995, 1997). Localised solutions after an oscillatory 
bifurcation have recently been examined in two dimensions by Riecke and Granzow 
(1998), who found that such solutions may be obtained in a pair of Ginzburg-Landau 
equations coupled to a mean-field equation. Of particular interest is their result 
that stable localised states can be found, as here, when the primary bifurcation is 
supercritical and the coefficients of the amplitude equations are real. 
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